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ABSTRACT 

RY Tau is a rapidly rotating Classical T Tauri star observed close to edge-on. The 
combination of new HST/STIS observations obtained in 2001 with HST/GHRS 
Archive data from 1993 has allowed us to get, for the first time, information on the 
thermal structure and the velocity law of the wind. The repeated observations of the 
Si III] and C III] lines show a lack of changes with time in the blue side of the profile 
(dominated by the wind contribution). Very high temperature plasma (logTg = 4.8) 
is detected at densities of 9.5 < logne(cm~^) < 10.2 associated with the wind. The 
emitting volumes are ~ (O.35-R0)^ suggesting a stellar origin. The wind kinematics 
derived from the profiles (Si III] , C III] and [O II] ) does not satisfy the theoretical 
predictions of MHD centrifugally driven disk winds. The profiles' asymmetry, large 
velocity dispersions and small variability as well as the small emitting volumes are 
best explained if the wind is produced by the contribution of several outfiows from 
atmospheric open field structures as those observed in the Sun. 
Keywords: stars: pre-main sequence — stars: winds, outfiows — stars : individ- 
ual(RY Tau) 



1 Introduction 

Winds from T Tauri Stars (TTSs) were assumed to be cold, in spite of the observa- 
tional evidence of our own Sun, because they were first detected as bipolar molecular 
outflows (Loren et al 1981) and optical jets with low ionization fractions. The detec- 
tion of C III] and Si III] semif or bidden line emission from RY Tau and RU Lup wind 
(Gomez de Castro & Verdugo, 2001; hereafter GdCV) provided the first evidence 
that electron temperatures as high as \ogTe{K) ~ 4.5 could be achieved at the base 
of the TTSs outfiows; this high temperature was derived from ultraviolet (UV) ob- 
servations and could not be inferred from the classical optical/infrared/radio tracers. 
Recently, Dupree et al (2005) have claimed that O VI absorption has been detected 
in TW Hya and T Tau. These new data represent a challenge to the current theories 
for outfiow generation in TTSs. In this letter, we provide observational evidence of 
hot winds being driven from the stellar atmosphere of RY Tau, a rapidly rotating 
TTS {vsini = 51.6 km/s, Hartmann & Stauffer, 1989). RY Tau has been classified 
as an UX Ori star: a pre-main sequence star which shows aperiodic eclipse-like min- 
ima caused by variable obscuration produced by circumstellar dust associated with 
the disk (Grinin, 1992). This suggests that the disk is seen close to edge-on. Also, 
the high projected rotation velocity of the stellar atmosphere suggests that the star 
inclination is close to 90". 

2 Data 

Some few spectral tracers have been selected for the subsequent study, namely, the 
C IV , C III] , [C II] , Si III] and [O II] lines observed in the UV spectrum of 
RY Tau (see Table 1 for more details on the observations and the line fluxes, and 
Fig. 1 for the proflles). All data are public and can be withdrawn from the HST 
Archive (they were obtained either by us or by other HST users). Data have been 
processed with the Routine Science Data Pipeline (RSDP). The major source of 
inaccuracy in the calibration is the centering of the target in the aperture that 
can account for as much as 20 km s~^ in the wavelength zeropoint accuracy of the 
GHRS data (sec GdCV) and <0. 1-0.2 pixels in the STIS data (e.g. velocity zeropoint 
accuracies better than 0.6 km s^^ and 2.8 km s^^ for the E230M and G230M spectra, 
respectively). WAVECALs were obtained during the observations and the peakup 
procedure (ACQ/PEAK) was used for the acquisition to center the source accurately, 
to about a 5% of the slit width (see www.stsci.edu/hst/stis/ for more details). This 
accuracy was cross-checked for the E230M spectra using as reference the narrow 
absorption features observed in the Fe ii, Mg ii and Mg i resonance lines produced 
by the interstellar medium. The STIS spectra were co-added after this check to 
improve the SNR of the data since no variations are detected. The profiles have been 
classified in Fig 1 by observing date. The left panel displays the very asymmetric 
C III] and Si III] profiles that led GdCV to conclude that these lines are produced at 
the base of the wind. The C IV line observed ~ 2 hours before is also displayed; the 
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profile is very similar to the C III] and Si III] lines. Bluewards shifted emission is 
also detected in the [O II] lines (see central panel) observed in Feb. 2001. The [C II] 
, Mg II and Mg I lines are obtained in the same spectra. A minimum is detected at 
a similar velocity in the Mg II lines; these lines have large optical depths in dense 
plasmas with T ~ lO^K. A narrow absorption component, probably of interstellar 
origin, is detected at the stellar rest velocity. There also seems to be an absorption 
component in the Mg I line at the same velocity. The last observations of the C III] 
and Si III] lines are displayed in the right panel of Fig.l. The profiles changed 
significantly from 1993 to 2001, although the bluewards shifted part of the profile 
remains stable. 



3 Interpretation and Analysis 

The detection of emission from high velocity gas in C IV , Si III] , C III] , [O II] , 
[O I] (Hamann 1994) indicates that there is a broad range of temperatures in the RY 
Tau wind at similar velocities. The wind properties seem to be rather stable; the 
emission peaks at similar bluewards shifted velocities in observations obtained over 
10 years, compare Fig. 1 with Hamann (1994) that measures a blueshifted [O I] 
emission peaking at -79 km/s. 

Line ratios can be used to constrain the wind temperature and density taking 
into account that the regions traced by each spectral tracer only overlap partially. 
Also occultation effects may be relevant (see Sect. 4). Line emissivities have been 
calculated for a collisional plasma over a broad range of densities {ue = 10'' — 
10^^ cm~^) and temperatures (logTe(i^) = 4.0 — 5.0) using the atomic parameters 
in the CHIANTI data base (Dere et al 1997) and assuming solar abundances and 
ionization equilibrium. 

The high temperature region: from the GHRS observations we derive: C IV /Si III] 
= 0.6 ±0.7 and Si III] /C III] =1.1 ±0.6. Then from Fig. 2 we derive: log(Te) ~ 4.8 
and 9.5 < log(?2e) < 10.2. Neither the line ratios nor the UV line fluxes vary by more 
than a factor of 2 in the fifteen observations available from 1979 to 2001 (GdCV). 
Notice that the high resolution data indicate that the variable contribution from 
0-velocity and redshifted material to the Si III] and C III] line fluxes does not 
significantly change their line ratios 0. 

The low temperature region: The [C II] emission feature consists of five spectral lines 
(see Fig. 1 and Table 1). The individual features are not resolved in RY Tau; a broad 
blend is observed instead. Some ratios among the individual lines in the feature are 
density sensitive in the range n^. = 10^ — 10^^ cm~^. An accurate determination of the 
electron density in the [C II] emission region requires us to derive simultaneously 
the ratios among the multiplet lines and the underlying [C II] profile. A least 
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squares scheme has been used to solve this inverse problem. The optimal fit gives 
a high lower limit of > 10^° cm~^ for the [C II] emitting region, similar to that 
determined from the C III] line. Note that the [C II] profile is too noisy to narrow 
this range (the contribution from the Fe II(uv3) line at vacuum wavelength 2328. llA 
is negligible as shown by the weakness of nearby features corresponding to the same 
multiplet). The fit also indicates that the underlying [C II] profile is similar to the 
C III] (STIS) profile; e.g. a blueshifted wind profile plus an excess from 0-velocity 
to the redwards shifted edge. Prom that fit, we also obtain ^Jfi 2326. 117) ^ 1.3 ±0.1. 
This ratio is temperature sensitive and from that we derive that 4.4 < log(Te) < 4.5. 

4 Implications for wind models 

TTSs outflows are characterized by their collimation at large scales, i.e.. it is ex- 
pected that whatever the ejection mechanism is, the wind terminal speed corre- 
sponds to a longitudinal motion along the collimation axis. Current models make 
use of magnetic fields to guarantee the large scale collimation. In this section we use 
the models of Gomez dc Castro & Ferro-Fontan 2005 (hereafter GdCFF) to compare 
to data. Disk wind models by other authors have a similar kinematics though the 
scales are different (see Uzdensky 2004 for a review). 

In GdCFF models, the wind is launched from the disk corona by a combination 
of thermal pressure and centrifugal force. This results in a radial expansion, away 
from the rotation axis, at the base and produces an increasing broadening of the 
profiles in this acceleration region. After this, the restoring toroidal component of 
the field narrows down the beam and the velocity field is progressively dominated 
by the motion along the rotation axis until it reaches the asymptotic regime, i.e. 
flow along the axis at the terminal velocity and basically no radial expansion. The 
evolution of the radial expansion, Vr, with the height above the disk is plotted in 
the top panel of Fig 3 for the warm disk winds model. Line emission is produced 
at different locations in the flow depending on densities and temperatures thus, the 
line profiles depend on the region of the wind sampled by each particular spectral 
indicator (as well as on the inclination). 

In order to compare these theoretical predictions with the observations, we have 
defined two observables that are measured only from the bluewards shifted emission 
of the relevant spectral tracers to avoid the pollution by the contribution of the 
accretion fiow. Notice that line emission from high temperature infalling plasma is 
produced at scales << i?* pumped by the X-ray radiation released at the accretion 
shocks. Thus occultation by the stellar disk causes the contribution to the line fiux 
from accretion to be redwards shifted (see Bcristain et al (2001) for more details). 
These observables arc the velocity dispersion, dV, and the line centroid,Vc. The line 
profiles may be understood as velocity-fiux plots {or V — f plots). The centroid of 

the bluewards shifted region of the profile is defined as, Vc = ^-^JnV^^^ and the 
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dispersion, dV, as, dV = w Jl^^ r for a uniformly sampled profile with 

V ^i=i 

step /SV and points. Vq represents the velocity at the blue edge of the profile and 
= ABSiVol /S.V) + 1. These observables have been calculated for three spectral 
tracers that sample very different plasma densities: Si III] , C III] and [O II] . The 
critical electron densities for the coUisional quenching of the C III] and [O II] lines 
are 1.5 x 10^° cm~^ and 0.5 x 10^ cm~^ respectively. For comparison, the critical 
density of the Si III] line is 2.0 x lO'^^ cm~^ for a fiducial temperature of 2.5 x lO^K. 
Vc and dV are plotted in the bottom panel of Fig 3 for the GdCFF theoretical 
predictions and various inclinations; notice that the [O II] line is predicted to be 
formed in the rarified expanding region while the GUI] and Si III] lines are predicted 
to be formed in the hot disk corona where the wind is accelerated. 

Vc and dV have been computed for the [O II] , C III] and Si III] profiles observed 
in RY Tau. To take into account the effect of the noise on these observables, we 
have computed these values both on the observed profiles and on filtered profiles 
obtained by applying a white noise filter to the data. As shown in Fig. 3, the major 
effect of the noise is to increase the velocity dispersion and shift the centroid some 
few km s~^ to the blue; however the overall trend remains. Two important results 
are derived from the figure: 

1. The line centroid shifts to higher velocity as the density traced by the spectral 
line becomes lower as expected (see critical densities), 

2. The velocity dispersion is high, > 40km s~^, and similar for all the tracers. 
In addition, the profiles are very asymmetric (see Fig. 1). Rather constant 
velocity dispersions and asymmetric profiles are expected for very small in- 
clinations, while for large inclinations the profiles should be symmetric and 
the dispersion increase steadily with the line centroid. However, all evidence 
points to RY Tau being observed does to edge-on (see Sect. 1). 

Thus RY Tau observations cannot be reconciled with models based on the centrifugal 
launching of the TTSs jets from the accretion disk. 

Another important piece of evidence comes from the emitting volumes. The 
volume of plasma producing the observed Si III] and C III] flux can be derived 
from the line emissivit}!^. The emitting volume producing the Si III] emission is 
2.07 X lO^^cm^ (z.e. (O.24R0)3) and similar to the derived for the C III] and C IV 
lines. Thus, line emission is produced in either a small structure or in dense clumps 
in the outflow. Notice that even adopting a filling factor of 0.1%-1% for the clumps, 
as that derived from observations of optical jets (Liseau et al 1996), the volume of 
the region producing the observed C III] and Si III] profiles just increases to about 
1.1 — 2.4 Rq. Thus, it is highly likely that mass ejection occurs from the stellar 
surface. 

^The emissivities are: j(C IV ) = 1.44 x lO^^erg cm"^ s"^ for log(re) = 4.8 and fog(ne) = H; 
j(Si III] ) = 1.49 X lO-^erg cm-^ for fog(Te) = 4.8 and fog(ne) 11; and j(C III] ) = 
0.97 X lO-^erg cm-^ s^^ for fog(Te) = 4.8 and \og{n^) = 11. 
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As RY Tau is a fast rotator the centrifugal lever arm from the stellar surface 
is an efficient source of thrust (see e.g. Sakurai 1987). Thus, the wind could be 
ejected from the stellar surface from open field structures like those observed in the 
Solar corona. The interaction between the wind and the circumstellar material (and 
a probable slow wind) could produce shocked knots like the Corotating Interac- 
tion Regions (CIRs) observed in the interplanetary medium and in young stars like 
AB Dor (Gomez de Castro 2002) that could contribute to the UX Ori phenomenon 
observed in RY Tau. 

This suggestive possibility would indicate that the observed profiles represent just 
an average of many small scale outflows and although these would not be stationary 
on a one-by-one basis, the average of all of them could produce a rather stable profile 
(provided that the spectral resolution is not good enough to resolve the individual 
components). If the magnetic field configuration is not axisymmetric (as observed 
in cool main sequence stars) this wind would naturally lead to the formation of 
nonaxisymmetric outflows. Tracking the variability of the C IV and C III] profiles 
is crucial for testing this possibility. 

Acknowledgments: We thank an anonymous referee for helping us to sharpen the 
definition of the observables in Fig. 3. 
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GHRS: DEC/31/1993 STIS: FEB/ZO/2001 STIS: MARCH/27/2001 
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Figure 1: UV lines observed in RY Tau and used in this work; RSDP processed 
data (see Table 1 for the binning) are plotted with a thin line and the 3-pixels 
average profile with a thick hne. The rest wavelength of the hnes and the velocity 
of the unresolved jet at ~ —80 km/s (from GdCV and Hamann 1994) arc marked 
with dashed lines. Left panel: Observations obtained in Dec. 31st, 1993 with the 
GHRS. The C IV line is blended with a very weak H2 feature (see Ardila et al 2002). 
Central panel: Profiles from STIS observations obtained in Feb. 20-21, 2001. The 5 
components of the [C II] multiplet are marked. Right panel: C III] and Si III] hnes 
observed in March 27th, 2001. Both lines show an excess of flux in the red wing 
compared with the 1993 observations; this excess is shaded in the figure. 
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Figure 2: C IV / Si III] versus Si III] / C III] flux ratios (in erg/s/cm-^) for 
collisional plasmas with electron densities between 10^ — lO^^cm^"^, temperatures 
log(Te)=4-5 and solar abundances. The value (with the error bars) derived from 
high resolution observations is marked with thick line; the 13 low resolution lUE 
observations are plotted with thin hnes. Isothermal and isochoric lines are plotted 
with continuous and dashed lines respectively. 
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Figure 3: Top: variation of the wind expansion velocity with respect to the jet axis, 
Vr, with the height above the disk, z, from GdCFF model The location of the 
Si III] , C III] and [O II] emission regions is sketched (see GdCFF for the detailed 
calculations). Bottom: Comparison between the theoretical predictions from the 
GdCFF model and RY Tau observations. Vc and dV values measured from the 
Si III] , C III] and [O II] profiles observed in RY Tau are marked with big squares, 
triangles and circles respectively. Vc and dV are plotted from the Si III] , C III] 
and [O II] profiles derived from the theory (small symbols) for various inclinations. 
High dV values are obtained from the observed profile while the slightly smaller dV 
values are measured over filtered profiles (see text). 
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Table 1: HST/STIS observations. 



Inst. 


Grating 


Res. 


Scaled"' 

A / clem 


Slit 


Spec. Range 

A 


Obs. Date 

(dd-mm-yyyy) 


Target^'') 
Features 


Flux '•"I 

lO^-*"** (crg/s/cm'^) 




G160M 


20000 


0.069 




1532 - 1567 


31-12-1993 


C IV 


F(C IV)= 8.2 ±7.1(9) 




G200M 


20000 


0.078 


LSA 


1880 - 1920 


31-12-1993 


Si III], C III] 


F(Si III])= 13 ±4 


STIS(9) 














F(C III])= 12 ±3 




E230M 


30000 


1/60,000 


0".2xO".2 


2303 - 3111 


19-02-2001 


[O II], [C II],Mg II, Mg I 


F([0 II])= 5.5 ± 1.4 




E230M 


30000 


1/60,000 


0".2x0".2 


2303 - 3111 


20-02-2001 


[O II], [C II], Mg II, Mg I 


F([C II]2326.l)= 19 ±2 




E230M 


30000 


1/60,000 


0".2x0".2 


2303 - 3111 


20-02-2001 


[O II], [C II], Mg II, Mgl 






G230M 


10500 


0.09 


52"x0".2 


1839 - 1929 


27-03-2001 


Si III], C III] 


F(Si III])= 26 ±3 
F(C III])= 15 ±3 



The element is a diode for GHRS spectra and a pixel for STIS spectra. 

The vacuum wavelengths of the target features arc: C IV(1548A, 1550A), Si III](1892A), C III](1908A), [O II] (247lA), 

[C II](2324.2A,2325.4A,2326.lA,2327.6A,2328.8A), Mg II (2796A,2803A),Mg I (2853A) 
(c) The Interstellar reddening law has been used since the extinction law towards T Tauri systems is not known. 

Parameters: Av=0.55mag; Aciv/^v = 2.60; Asmi/Av = 2.66; Acm/Av = 2.66; Acii/Ay = 2.60; Aon/Ay = 2.44. 

The flux of the [C II] (2326. lA) has been derived from the fitting of the [C II] multiplet (see text). 
^'^^ Goddard High Resolution Spectrograph removed in February 1997 
('^^ Large Science Aperture (2") 

The H2 (3-8 R(3)) line contribution to the C IV flux is negligible (see Ardila et al 2002). 

Space Telescope Imaging Spectrograph stopped science operations in August 2004 
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